void draw(){
  TFile *f = new TFile(getenv("TMP1"), "UPDATE");
  TTree *etf = (TTree *)f->Get("etf");
  TTree *et = (TTree *)f->Get("et");
  etf->AddFriend(et);
  char cmd[512]{};

   TCutG *cutg = new TCutG("CUTG",16);
   cutg->SetVarX("aoz0+0.00099*pdcTaU[1][0]-0.0031*x_f1");
   cutg->SetVarY("z0");
   cutg->SetTitle("Graph");
   cutg->SetFillStyle(1000);
   cutg->SetPoint(0,2.11938,36.5816);
   cutg->SetPoint(1,2.09012,36.9224);
   cutg->SetPoint(2,2.08051,36.5065);
   cutg->SetPoint(3,2.08639,35.981);
   cutg->SetPoint(4,2.10805,35.6922);
   cutg->SetPoint(5,2.12296,35.5074);
   cutg->SetPoint(6,2.13358,35.6865);
   cutg->SetPoint(7,2.14276,35.7789);
   cutg->SetPoint(8,2.16384,35.4612);
   cutg->SetPoint(9,2.17933,35.5248);
   cutg->SetPoint(10,2.18908,35.8828);
   cutg->SetPoint(11,2.18607,36.2986);
   cutg->SetPoint(12,2.16126,36.5527);
   cutg->SetPoint(13,2.14175,36.7491);
   cutg->SetPoint(14,2.12684,36.6047);
   cutg->SetPoint(15,2.11938,36.5816);

	TCanvas *c = new TCanvas("c", "c", 1200, 800);
	c->SetGrid();
	c->ToggleToolBar();
	c->ToggleEventStatus();
  for(int i = 0; i < 10; i++){
    double kf1 = 0.001 + i*0.0001 - 0.0005; // 0.001-0.002
    for(int j = 0; j < 10; j++){
      double ketf = 0.003 + i*0.0001 - 0.0015; // 0.001-0.002
      sprintf(cmd, "aoz0+%f*pdcTaU[1][0]-%f*x_f1>>h(1000, 2., 2.2)", ketf, kf1);
      etf->Draw(cmd, "pdcTaU[1][0]!=-999.&&x_f1!=-999.&&x_f1!=9&&x_f1!=25&&CUTG");
      c->Update(); c->Modified();
      getchar();
      // sleep(1);
    } // end for over j
  } // end for over i
} // end main function